Possible mechanism for glass-like thermal conductivities in crystals with off-center 

atoms 
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In the filled Ga/Ge clathrate, Eu and Sr are off-center in site 2 but Ba is on-center. All three 
filler atoms (Ba,Eu,Sr) have low temperature Einstein modes; yet only for the Eu and Sr systems is 
there a large dip in the thermal conductivity, attributed to the Einstein modes. No dip is observed 
for Ba. Here we argue that it is the off-center displacement that is crucial for understanding this 
unexplained difference in behavior. It enhances the coupling between the "rattler" motion and the 
lattice phonons for the Eu and Sr systems, and turns on/off another scattering mechanism (for IK < 
T < 20K) produced by the presence/absence of off-center sites. The random occupation of different 
off-center sites produces a high density of symmetry-breaking defects which scatters phonons. It 
may also be important for improving our understanding of other glassy systems. 

PACS numbers: 



An important, universal, characteristic of glass-like 
systems is a low thermal conductivity k; k(T) varies 
roughly as for T < IK, forms a plateau/dip region 
somewhere in the range 5-40K, and then increases slowly 
at higher temperatures. The low T behavior (< IK) is 
well understood, and is attributed to a broad distribu- 
tion of tunneling centers (or two level systems), both 
in glasses and in some disordered crystals that exhibit 
glass-like behavior such as (KBr)i_a;:(KCN)afi, while at 
high temperatures, Rayleigh scattering becomes impor- 
tant. However the intermediate plateau/dip region is less 
well understood. Several explanations for the plateau 
have been proposed, but the model of Grannan etal^ for 
the mixed KBriKCN system has an important general 
mechanism and appears to be applicable for both glassy 
and disordered crystalline materials. In their model, the 
plateau/dip is produced by a nearly localized mode [the 
libration modes of the CN ion about its center of mass, at 
THz frequencies] that resonantly scatters phonons very 
effectively. Similar models have been applied to glassy 
selenium:^ and disordered garnets^. The same approach 
has also been used to describe the glass-like thermal con- 
ductivity in the off-center Eu and Sr filled clathratea^; in 
this case a "rattler" mode with a low Einstein tempera- 
ture in the range 30-150K (i.e. at THz frequencies) plays 
the role of the nearly localized mode. However this model 
fails in explaining the difference in behavior between the 
Ba and Eu/Sr filled clathrates - no dip is observed for the 
on-center Ba system. It also cannot explain the k cx T 
dependence observed from roughly I-IOK^. Here we will 
argue that it is the off-center displacement of Eu and Sr 
that leads to the plateau/dip region in these systems. 

In several compounds with large unit cells (skutteru- 
dites and type I clathrates), large cages or voids oc- 
cur in the structure which can be "filled" with several 
types of atoms. This dramatically alters their physical 
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When the "filler" ion is considerably 



smaller than the void, it is loosely bound and can "rat- 
tle" around - hence its name. The weak binding leads to 
a nearly local modeiS described by a low Einstein tem- 



perature. A low thermal conductivity has been found in 
a number of such systems^i^*ii and a glass-like model for 
k(T), used for the Eu and Sr clathrates^. Note that a 
very low value of n is crucial for thermoelectric appli- 
cations because the figure of merit ZT = TS^cte/k (S is 
the Seebeck coefficient, ae, the electrical conductivity) 
can be significantly increased. Hence understanding the 
mechanisms that lead to a small k, and particularly their 
ranges of validity, is very important. 

Here we consider k{T) for the (type I) Ga/Ge 
clathrates whose structure is formed of two cages with 
internal voids large enough to house another atom (sites 
1 and 2); atoms placed in them, such as Eu, are re- 
ferred to as Eul and Eu2. In the larger site 2 cage, Eu2, 
Sr2, and Ba2 all form rattlers with quite low Einstein 
temperatures^; however Ba2 is on-center while both Eu2 
and Sr2 move off-centeri2iiiiiiiS*iSJI. 

In Fig. we re-plot the thermal conductivity data 
of Nolas eta&; k is small for Eu and Sr, and there is a 
well defined dip associated with the Einstein tempera- 
tures of the rattlers'*. Nolas etal model the low thermal 
conductivity below 20K for the Eu and Sr compounds 
using a broad distribution of off-center tunneling states. 
At slightly higher temperatures, they used two resonant 
scattering terms (two Einstein temperatures) to explain 
the dip (near 25K), plus a Rayleigh scattering term which 
dominates at high T. Using this sum of terms, they were 
able to fit k(T) for both EugGaigGeao and SrgGaieGeso 
over the range 5- lOOK. However a problem arises in us- 
ing a tunneling model with tunneling splittings up to ^ 
15K;- first the inferred distribution of tunneling states 
from ultrasonic measurements is much narrower than in 
a glassi^, and second (more importantly), very recent 
Mossbauer measurements indicate that the tunneling fre- 
quency is at most 0.44 Ghz - i.e. ~ 0.02K*^. If one cuts 
off the tunneling density of states with a step function 
even at IK, then k rises too rapidly below 15K. 

Similar experimental results for k{T) of these filled 
clathrates have also been observed by Sales eiafe (See 
Fig. Qd); however, at low T, k is considerably lower and 
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FIG. 1: a) The lattice thermal conductivity for EugGaigGeso 
and SrgGaieGeao from ref 0- Solid lines - the new fits 
discussed below; dashed line - the calculated result for 
SrgGaieGeso if the resonant scattering contribution were re- 
duced by a factor of 10. It shows a dependence similar to 
that of BagGaieGeso in b. The inset shows a log-log plot, b) 
Similar data for the Eu Sr, and Ba clathrates from ref|3. Solid 
lines show the fits; for Eu and Sr the same model is used as 
for a) but G and the d are larger; for Ba, the d are small 
and a finite sample size and an Umklapp term are needed at 
low T. 



the dip is less pronounced. More importantly, for the Eu 
and Sr data, k(T) decreases nearly linearly with T below 
lOK, which can't be fit with the above model. 

The T-dependence for BagGaigGeao in which Ba2 is 
on-centen^, is more perplexing (See Fig. ^p). For this 
system, n is much higher overall and there is no clear 
dip in k(T), although the Einstein temperatures arc only 
slightly larger than those of Eu and Sr^; thus the rattler- 
phonon coupling must be greatly reduced^". Also, the 
T-dependence below 10-12K is much faster (k ~ T'^) for 
the Ba clathrate. These results raise two important issues 
- what determines the strength of the coupling between 
the rattlers and the phonons and how can one understand 
the T-dependence of k from 1-1 5K? We propose another 
mechanism to explain this temperature range; it may also 
contribute to k in other glassy systems. 

The off-center displacement appears to play two ad- 
ditional, crucial roles in increasing the phonon scatter- 
ing (in addition to producing tunneling centers). First 



from XAFS, the off-center atom is partially bonded to 
the nearest neighbor Ga/Ge atoms on the side of the 
cageii, forming a random array of symmetry-breaking 
defects. We show below that a high density of such de- 
fects can lead to a k ^ T dependence at low and medium 
temperatures. A likely more important consequence of 
this bonding is that the motions of the off-center rat- 
tler atom are coupled more strongly to just a few of the 
closest atoms in the surrounding Ga/Ge framework and 
hence to the phonon modes - i.e. the motions of the Eu2 
and nearest Ga/Ge atoms are positively correlated. 

We first address the use of the tunneling model and 
the range of reasonable tunneling splittings. Sethna and 
coworkers^ suggest that for the light CN molecule the 
tunneling splittings (for reorientations of the CN axis) 
might extend to about IK; measured tunneling split- 
tings on similar, more dilute systems (CN~ and 0H~) 
are less than SK^i*^. For heavy off-center ions such as 
Ag+ and Cu+ in the alkali halides, the observed tunneling 
energies are also low, much less than IK^S; only for the 
light off-center ^Li+ ion does the tunneling splitting reach 
1.2K at normal pressures^'^. The tunneling behavior can 
be finely tuned by applying hydrostatic pressure to re- 
duce the potential barrier between equivalent off-center 
sitesSiS^, but even in such cases the tunneling splittings 
for heavy atoms such as Ag+ are about IK when the off- 
to on-center transition begins to occur'^'*. 

Within the simple two-well tunneling model, having 
a tunneHng splitting (^ 10-15K) that is ~ 10% of the 
attempt frequency (the rattler frequency) is inconsistent 
with the assumptions for the model; for a heavy mass 
such as Eu, it leads to an unreasonably low potential 
barrier between the wells (< 1 meV (12K); - i.e. less 
than the tunneling splitting). On the experimental side, 
a broad distribution of tunneling states should also lead 
to a linear contribution to the heat capacity at low T 
when the tunneling states are excited; however the heat 
capacity for the Ba, Sr, and Eu systems at low T is expo- 
nential and can be fit to an Einstein model^ ; thus there is 
no evidence for tunneling states extending up to 10-15K 
in such data. With the recent observations that the tun- 
neling frequencies are at most 0.02KiS for the Eu system, 
another mechanism is needed to describe k at low T. 

Since atomic forces are generally a very strong func- 
tion of the distance between the atoms/ions, the force 
between the central atom (e.g. Sr2 or Eu2) and each of 
the nearest equivalent Ga/Ge atoms in the cage, will in- 
crease substantially for an off-center displacement. The 
resulting vibrations of the off-center rattler/cage system 
then depend on the rigidity of the cage (the Ga/Ge cage 
is quite stiff^'^^), the Eu2-Ga/Ge bond strength, and the 
number of Eu2-Ga/Ge bonds. In the limit of an idealized 
rigid cage, there would be no motion of the cage atoms 
(i.e. no phonon-coupling) as the rattler atom vibrates; 
the reduced mass of the local mode would be that of the 
rattler atom as assumed previousljii^*iSiii. For the oppo- 
site extreme of a very soft cage, the local mode reduced 
mass for motion along the bond direction would approach 
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the reduced mass of the atom-pair involved, e.g. Eu and 
Ge in EugGaigGeso. In this case, there would be large 
motions of the nearest cage atoms (Ge) as a result of 
a local mode vibration, which would couple directly to 
the phonons of the clathrate framework. Because the 
clathrate cage is quite stiff, the reduced mass is expected 
to be close (but not equal) to the free rattler mass. 

To further understand this coupling, consider a 1-D 
model of a cage (mass M) connected to a rattler (mass 
m), and let the respective displacements be x and X. 
Then Mx = mX; we expect the matrix element for scat- 
tering to be proportional to x/X=m/M, and the cou- 
pling to {x/X)^. The effective mass of the cage decreases 
when the rattler moves off-center (because only a small 
fraction of the cage atoms are directly bonded to the rat- 
tler) and thus the coupling to phonons increases. Also 
note that for the off-center case, rattler vibrations could 
be either ~ radial (along the bond) or perpendicular to 
the bond. The latter are similar to libration modes and 
would have a lower Q e than for radial motion which was 
not considered for the KChKCN system. 

To investigate the various contributions to n we first 
reproduced the calculation of Nolas et alf^-. k is given by: 
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1{lo) = (A^/ + A-,1 + A^;^)-! +;™„ 
D{hLu/k)^ 



Aihuj/kB )tanh{huj /2kBT) 
+ {A/2){kB/hLj + B-^T-"^)-^ 
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where Cg is the speed of sound, C(a;,T) is the heat ca- 
pacity of phonons with frequency to and l(uj) is the total 
mean free path which has three components - Rayleigh 
scattering (A^j), resonant scattering from the Einstein 
modes (A^es) and a contribution from a broad distribu- 
tion of tunneling states (A^^); the lower limit is con- 
strained to lmi7i- The Einstein frequencies were ob- 
tained from structural data and the constants A, B, Ci, 
C2, D, 7_Bi and Imin are given in ref. 0. We obtained 
similar fits although the constants varied slightly. We 
then reduced the coupling to the resonant modes - the 
Ci constants - by a factor of 10 and left the other terms 
unchanged. This result is shown by the dashed line (Fig. 
^i) and is very similar in shape to the data of Sales^ for 
BasGaieGcao (Fig^jD). Thus a major part of the dif- 
ference in the T dependence of k between BasGaieGcao, 
and the Eu or Sr compounds can indeed be explained by 
a decreased coupling between on-center Ba2 and phonons 
in the Ga/Ge framework as described above. 

An important parameter for probing the rattler- 
phonon coupling is the reduced mass fi, of the local mode 
oscillator: however, fi is not an easy parameter to obtain. 
One recent study of rattlers in the skutteruditesS- , found 



the reduced mass to be somewhat less than the atomic 
mass of the rattler (error ~ 20%), but more measure- 
ments are needed. 

A second consequence of the rattler being partially 
bonded to a side of the cage is the formation of a high- 
density, random network of symmetry-breaking mass de- 
fects - i.e. one of the four off-center orientations is occu- 
pied, but the off-center site varies randomly throughout 
the crystal. The disorder introduced by the random oc- 
cupation of one of the off-center minima provides insight 
as to why the crystal exhibits glass-like behavior. This 
model has similar features to the disorder introduced via 
irradiation of a quartz crystal^LS but there is a key dif- 
ference - the atoms are randomly displaced in the irra- 
diated sample which leads to a range of bond lengths 
and hence to the possibility of a broad distribution of 
tunneling states. For the Eu/Sr clathrates the off-center 
displacement is well defined; the disorder arises from a 
random occupation of the four off-center sites. Locally 
such disorder changes the speed of sound as the phonon 
wave passes a rattler atom; the resulting phase changes 
must then be included. Note that in this case there is 
an interference in the scattering between different rat- 
tler sites, whereas Rayleigh scattering assumes a single 
scattering process. 

Ziman has considered this type of scattering in some 
detaiP^ - in his model, the effective mean free path for 
this contribution is 
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where q is the wavenumber, Scg is the variation of the 
sound speed, L is the smallest distance between scatter- 
ers, and ^ is a numerical constant of order 1. For very 
short wavelengths A << L, ^disorder goes to a constant - 
but that is already included in the model of Nolas et alA 
by the Imin term. Using the Debye model (lu ~ Cgq) we 
then replaced the tunneling contribution in Equ. |21 
by 



^disorder - Guj"^ . 



(5) 



In Fig. HJi we show that with Geu ^ 3.7 10"^^s^/m 
{Gsr = 6.4 10~^°s^/m), and a shght change of the Ci 
parameters this model describes the T behavior of k (of 
Ref. 0) just as well as by using the tunneling model. 
Taking L ~ 5 10~^'^m, these numbers imply Scg/cs ~ 
7% (Eu) and 3%(Sr). Note that this disorder mechanism 
is important at moderately long wavelengths - i.e. for 
temperatures up to 20-30K, and alone would yield k(T)^ 
T. 

Attempts to fit the Sales et al data for Eu and Sr— using 
the model of Nolas etai^ were not successful^ - the data 
from 2-lOK vary nearly linearly with T, and cannot be 
fit by the ~ T^ dependence for the tunneling model at 
low T. However, we can fit the data quite well (Fig. ^p) 
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using A^-^Q^^g^ but larger values for G and the d are 
needed to fit the lower thermal conductivity. 

For Ba, very much smaller values for Ci are needed, 
which confirms the earlier assertion that the coupling be- 
tween phonons and the rattler vibrations is reduced. In 
addition, the ^^isorder term cannot fit the low T data for 
Ba below 12K; k is higher and the approximately de- 
pendence suggests that a finite sample size contribution 
is needed at low T - details will be given in a separate 
paper—. The lack of a A^.^^^^^^ term is consistent with 
Ba being on-center, as in that case the A^-^^^^^^ term 
should disappear. 

In summary, we have pointed out that an off-center 
rattler enhances the scattering of phonons in two ways. 
First, the off-center displacement increases the coupling 
of the rattler motion to the phonons and hence increases 
a dip in the thermal conductivity near ~30K. This pro- 
vides an explanation for the the lack of a dip in k(T) 
for the on-center Ba system. Second, the off-center atom 
on the side of the cage, introduces a quasi-random set 
of "extra" atoms that will produce local changes in the 
sound velocity. This in turn leads to a mean free path 
that varies as w"^, which provides an alternative expla- 
nation for the reduction of the thermal conductivity for 
IK < T < 20K. It results in a linear T dependence in 



this range which fits the data of Sales etaP well. 

Several questions remain to be answered - for the 
clathrates, what are the effective masses of the local Ein- 
stein modes and what are the actual tunneling splittings 
(there will be 3 levels) in these systems? Can the Sr sys- 
tem be moved on-center via high hydrostatic pressure and 
thus allow an investigation of variations in the rattler- 
phonon coupling? Can these systems be optimized for 
thermoelectric applications over a particular T-range, by 
varying the off-center rattler's mass and/or Einstein tem- 
peratures? More generally does the off-center displace- 
ment play a role in producing the plateau/dip regions in 
other glassy systems? 
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